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Abstract 

This  investigation  considers  the  mode -conversion  process 
encountered  in  the  lower-hybrid  heating  scheme  at  the  edge  of  the 
plasma,  where  the  frequency  of  the  external  radiation  is  compar- 
able to  the  local  value  of  the  electron  plasma  frequency.  The 
problem  is  formulated  in  terms  of  a modulational  representation 
of  the  electric  field  that  permits  the  description  of  the  full 
space-time  evolution  of  the  process.  The  numerical  study  of  the 
model  equations  confirm  that  a backward  wave  is  excited  at  the 
edge  and  that  it  leads  to  a net  energy  flow  directed  toward  the 
interior  of  the  plasma.  For  small  external  power  levels,  the 
flow  is  characterized  by  the  propagation  of  a leading-edge  pulse 
whose  arrival  at  a given  spatial  location  marks  the  onset  of  the 
steady-state.  At  large  external  power  levels,  the  self-consistent 
modification  of  the  density  profile  by  the  ponderomotive  force  is 
found  to  quench  the  mode-conversion  process,  thus  causing  a signifi- 
cant reduction  in  the  amount  of  RF  energy  that  can  be  coupled  from 
the  external  source  to  the  interior  of  the  plasma. 
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Abstract 

This  investigation  considers  the  mode-conversion  process 
encountered  in  the  lower-hybrid  heating  scheme  at  the  edge  of  the 
plasma,  where  the  frequency  of  the  external  radiation  is  compar- 
able to  the  local  value  of  the  electron  plasma  frequency.  The 
problem  is  formulated  in  terms  of  a modulational  representation 
of  the  electric  field  that  permits  the  description  of  the  full 
space-time  evolution  of  the  process.  The  numerical  study  of  the 
model  equations  confirm  that  a backward  wave  is  excited  at  the 
edge  and  that  it  leads  to  a net  energy  flow  directed  toward  the 
interior  of  the  plasma.  For  small  external  power  levels,  the 
flow  is  characterized  by  the  propagation  of  a leading-edge  pulse 
whose  arrival  at  a given  spatial  location  marks  the  onset  of  the 
steady-state.  At  large  external  power  levels,  the  self-consistent 
modification  of  the  density  profile  by  the  ponderomotive  force  is 
found  to  quench  the  mode -conversion  process,  thus  causing  a signifi- 
cant reduction  in  the  amount  of  RF  energy  that  can  be  coupled  from 
the  external  source  to  the  interior  of  the  plasma. 
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I.  Introduction 

The  goal  of  an  auxiliary  plasma  heating  scheme  aimed  at  achieving 
thermonuclear  temperatures  ( > 10  KeV)  consists  of  delivering  a large 
amount  of  external  energy  to  the  interior  of  a pre-existing  plasma  having 
a moderate  initial  temperature  (.  <3  KeV).  The  various  auxiliary  heating 
schemes  being  considered  at  the  present  time  differ  mainly  in  the  form 
of  the  primary  energy  used,  e.g.,  the  neutral-beam  heating  approach  uses 
the  streaming  or  kinetic  energy  of  the  injected  particles,  while  the  RF 
heating  schemes  utilize  electromagnetic  energy  directly. 

Regardless  of  the  particular  heating  method  considered,  in  its  im- 
plementation phase  one  always  encounters  the  following  major  physics 
questions:  1)  How  is  the  external  energy  transferred  to  the  ions  once 

it  reaches  the  interior  of  the  plasma?,  2)  How  is  the  external  energy 
transported  from  the  edge  of  the  plasma  to  the  interior?,  and  3)  How 
does  the  external  source  couple  to  the  surface  of  the  plasma?  It  should 
be  noted  that  the  existence  of  potentially  deleterious  answers  to  any  of 
these  questions  may  render  the  entire  heating  scheme  inapplicable  as  far 
as  the  attainment  of  fusion  conditions  is  concerned,  even  though  the  scheme 
may  be  quite  meritorious  on  the  grounds  of  the  existence  of  a well-devel- 
oped technology  or  because  of  its  ease  of  implementation. 

Historically,  it  has  been  the  case  that  the  interest  of  the  plasma 
physics  community  has  been  focused  mainly  on  the  issues  surrounding  ques- 
tion 1)  and  only  peripherally  on  those  matters  pertaining  to  questions  2) 
and  3).  However,  with  the  ever  increasing  need  to  assess  the  future  poten- 
tial of  auxiliary  heating  schemes  it  becomes  important  to  consider  the  funda- 
mental physics  of  processes  that  may  be  encountered  in  answering  questions 


I 
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2)  and  3).  It  is  in  this  general  direction  that  the  present  investigation 
is  aimed.  Specifically,  this  study  isolates  some  of  the  surface-coupling 
properties  of  external  electromagnetic  radiation  in  the  lower-hybrid  fre- 
quency range. 


The  particular  case  which  pertains  to  the  present  investigation  is  the 

lower-hybrid  heating  scheme  proposed  by  Stix1  over  a decade  ago.  A large 

fraction  of  the  theoretical  effort  in  this  topic  has  been  devoted  to  elluci- 

dating  the  numerous  processes  which  may  occur  deep  in  the  interior  of  the 

plasma,  where  the  frequency  of  the  external  radiation,  go,  matches  the  peak 

value  of  the  lower-hybrid  frequency  gotu  = w . / [1  + (go  /C2  )2]'2,  in  which  io 

LH  pi  1 pe  e 1 pe 

and  ui  ^ refer  to  the  electron  and  ion  plasma  frequencies,  and  represents 
the  electron  gyrof requency . 


Over  the  past  year,  it  has  become  increasingly  evident  that  major 
difficulties  may  be  encountered  in  this  heating  scheme  in  the  process  of 
transporting  the  RF  energy  from  the  edge  of  the  plasma  to  the  interior, 
i.e.,  before  the  condition  go  = go  is  ever  attained.  Numerous  theoretical 

Ln 

2-4 

studies  of  parametric  instabilities  , of  lower-hybrid  cone  filamentation 

7 

as  well  as  explicit  experimental  measurements  of  filamentation  phenomena 

and  associated  plasma  heating  in  small-scale  experiments,  have  contributed 

to  the  general  realization  of  this  point  of  view.  Yet,  a more  direct  and 

compelling  argument  has  been  provided  by  the  results  of  the  large-scale 

8 

lower-hybrid  heating  experiments  performed  in  Alcator  at  MIT,  and  in  the 

9 

ATC  tokamak  at  Princeton  . 

Although  the  details  of  the  outcome  of  these  experiments  are  difficult 


< 


to  interpret,  it  is  important  to  recognize  that  both  of  these  studies  are 
consistent  with  the  view  that  the  external  RF  energy  is  efficiently  injected 
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from  the  source  waveguide  into  the  plasma-filled  vacuum  vessel,  and  there- 
after it  propagates  in  the  vicinity  of  the  plasma  edge,  with  only  a rela- 
tively small  fraction  of  the  injected  energy  penetrating  to  the  interior. 
Therefore,  it  is  of  interest  to  isolate  some  of  the  basic  processes  which 
may  play  a role  in  the  coupling  of  external  RF  energy  in  the  lower-hybrid 
frequency  range  at  the  edge  of  a plasma. 

The  present  study  is  concerned  explicitly  with  the  mode-conversion 
process  which  is  encountered  in  the  low  density  region  at  the  edge  of  the 
plasma  where  the  frequency  of  the  incoming  radiation  is  of  the  order  of  the 

local  electron  plasma  frequency,  i.e.,  u>  - u .To  simplify  the  problem 

pe 

a slab  model  is  considered  in  which  the  density  varies  in  the  direction 

perpendicular  to  the  external  magnetic  field.  Unlike  previous  studies  of 

10,11 

the  mode -conversion  process  , which  have  dealt  exclusively  with  the 
steady-state  properties,  the  present  formulation  uses  a modulational  repre- 
sentation for  the  electric  field  which  permits  the  description  of  the  full 
space-time  evolution  of  the  detailed  process. 

The  numerical  study  of  the  relevant  model  equations  of  the  problem 

demonstrate  that  in  the  region  where  w > u a backward  wave  is  excited 

pe 

and  leads  to  the  propagation  of  RF  energy  away  from  the  plasma  edge  toward 
the  interior.  An  interesting  feature  of  this  transport  is  the  existence  of 
a propagating  leading-edge  pulse  whose  arrival  at  a given  spatial  location 


marks  the  onset  of  the  steady-state  pattern. 

The  present  formulation  also  permits  the  study  of  more  complicated  non- 
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linear  processes  which  may  take  place  at  the  plasma  edge  when  large  ex- 
ternal power  levels  are  applied.  In  particular,  the  present  study  has 
considered  the  effects  produced  by  the  profile  modification  caused  by 
the  ponderomotive  force  associated  with  the  intense  electric  fields.  It 
is  found  that  this  effect  can  alter  the  mode -conversion  process  to  a large 
degree,  thus  leading  to  significant  reduction  in  the  penetration  of  the 
external  RF  energy  into  the  plasma. 

The  manuscript  is  organized  as  follows:  Sec.  II  presents  the  formulation 

of  the  problem,  its  properties,  and  associated  conservation  laws.  Sec.  Ill 
discusses  briefly  the  boundary  conditions  and  numerical  scheme  used  to  study 
the  model  equations.  Sec.  IV  presents  the  results  of  the  time  evolution 
of  mode -conversion  in  both  the  linear  and  nonlinear  regimes.  A summary  of 
the  results  obtained  is  presented  in  Sec.  V. 
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II.  Formulation 

To  simplify  the  geometry  of  the  present  problem  we  consider  ex- 
plicitly a plasma  slab  whose  zero-order  density  nQ(x)  increases  along 
the  x-direction,  which  is  taken  to  be  perpendicular  to  the  external 
magnetic  field  Bq,  as  is  illustrated  in  Fig.  1.  The  z-axis  of  the  carte- 
ssian  coordinate  system  points  along  Bq,  while  the  y-axis  is  perpendicular 
to  the  x-z  plane,  and  constitutes  the  ignorable  coordinate  of  the  problem, 
i.e.,  there  are  no  explicit  functional  dependences  on  it.  The  total 

electric  field  £ is  expressed  in  the  form 
' 

£(x,z,t)  = ECx,t)ei(-^z  ' + c.c.  (1) 

In  Eq.  Cl)  k represents  the  wavenumber  along  the  z coordinate  and 
is  assumed  to  be  fixed,  physically,  by  a suitably  designed  external  slow- 
wave  structure  or  simulated  by  a properly-phased  waveguide  array  Cgrill). 
The  fundamental  role  played  by  such  external  devices  is  that  they  permit 
the  experimenter  to  choose  the  proper  value  of  k that  optimizes  the 
penetration  of  radiation  into  the  plasma.  In  particular,  these  devices 
permit  the  generation  of  radiation  such  that  k > kQ  = oj/c,  where  c is  the 
speed  of  light,  and  thus  kQ  represents  the  vacuum  wavenumber.  As  is 
well  known1 ’ 10,11  this  prerequisite  must  be  met  by  the  external  radiation 
in  order  to  satisfy  the  lower-hybrid  accessibility  criterion. 

In  the  context  of  the  present  problem  w is  equal  or  comparable  to 
the  value  of  the  peak  lower-hybrid  frequency  of  the  plasma;  such  a value 
is,  of  course,  attained  deep  in  its  interior.  In  contrast,  the  region  of 
interest  to  the  present  calculation  is  that  where  co  is  comparable  to  or 
within  an  order  of  magnitude  smaller  than  the  local  value  of  the  electron 
plasma  frequency,  i.e., 
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U)  = (UtU)  , 'V  - ) , 

LH  peak  pi  peak 


thus,  the  region  of  interest  is  defined,  for  a given  density  profile,  by 


$ «>  (x)  % 10u 

rc 


Maxwell's  equations  for  the  radiation  fields  £ and  B are 


O D 4TT  . 1 3 P 

'",  = 7'*;se 


which  can  be  combined  to  yield 


_ _ g 4i  3 . 

V x V *£  = --j.  — J 


1 3 2 ? 


considering  Eq.  (1)  and  making  use  of  the  modulational  approximation 
(i.e.,  neglecting  the  32E/3t2  term)  reduces  Eq.  (5)  to 


V x v x E = 


k2 

4 it  iut  3 . ,2C  -■  0 3 _ 

7 e 7T  J ~ + 2l  — rr  E 

ot-  O ^ /*»  nr 


ui  3t  _ 


In  the  spirit  of  the  modulational  representation  we  take  the  self-consistent 
RF  current  in  the  plasma,  j_,  to  be  determined  by  the  instantaneous  value 
of  the  electric  field,  i.e.,  we  neglect  time-history  contributions  to  the 
linear  dielectrics.  Accordingly, 


4tt  iut  3 . r 

—*■  e r—  j + E = 

it)*-  3t 
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where,  e is  the  cold  plasma  dielectric  tensor  having  the  form 


e e 

xx  xy 

£ £ 

yx  yy 


0 0 £ 


0 

0 

zz 


(8) 


The  ordering  of  frequencies  in  the  lower-hybrid  scheme  is  such  that 
<<  id  « thus  the  components  of  V"  are 

2 


0) 
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-EL 

yy 
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xy 
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(9) 


Using  Eq.  (7)  in  Eq.  (6)  yields 


— |-  E + V * E = r-i-  V x v x E 
U)  3t  kz 

0 


(10) 


Writing  out  explicity  the  V x V operator  in  Eq.  (10)  produces  the  following 
coupled  partial  differential  equations  for  the  components  of  the  electric 
field 


2i  3 . - 1 3 ^ c r /•  k ^ 2 1 c c 

“ 3t  Ey  + kf  3^2  Ey  + Ieyy  “ ("1T  ] / = ' VEx 


(11) 


2i  3 „ 1 32  _ ^ c . k 

a)  3t  z + k 3x2  z + Ezz  z _ 1 k 2 
o 0 


— E 

3x  x 


(12) 


2i  3 , r 
oi  3t  Ex  + ^exx 


] E = i r-r  E + £ E 
1 x k-  3x  z v-"  v 
o 


XV  X 


(13) 
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0 


(17) 
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Eq.  (17)  has  the  general  form  of  a damped  harmonic  oscillator  equation 
with  variable  coefficients..  The  contribution  from  the  term  proportional 
to  gives  rise  to  a damping  factor  which  is  important  only  for  ex- 

tremely sharp  gradients  in  the  dielectrics.  The  coefficients  multiplying 
the  term  proportional  to  are  of  crucial  importance  to  the  behavior  of 
the  system  because  their  combined  sign  determines  whether  or  not  the 
radiation  can  penetrate  into  the  plasma.  In  those  regions  of  extremely 
low  density  where  wpe00  < u>  it  follows  that  z < 0.  However,  as  one 
moves  across  the  = w point,  e^  = 0;  finally,  as  one  proceeds  to 

the  interior  of  the  plasma  where  upeW  > w,  z^z  > 0.  Accordingly,  the 

determining  factor  for  the  radiation  to  propagate  beyond  the  w = u (x) 

pe 

point  is  that  e/e^  < 0.  Since  in  this  region  e = 1 this  implies  that 
c <0,  or 


X2  = (p-)2  > 1 


If  this  condition  is  satisfied  then  the  external  radiation  mode-converts 
into  an  internal  plasma  oscillation  which  proceeds  to  propagate  to  the 
interior.  However,  the  mode-converted  oscillation  can  still  be  reflected 
at  the  point  where  e = 0.  Hence,  if  the  mode-converted  wave  is  ever  to 
reach  the  lower-hybrid  resonance  this  implies  that 


1 U) 

»2  ■ (|-)2  > i * 


for  all  the  density  values  up  to  the  resonance  layer.  Eq.  (19)  is  what 

one  refers  to  as  the  accessibility  criterion,  which  for  the  typical 
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parameters  of  relevance  to  thermonuclear  fusion  implies  that  (k/kQ)2  > 2, 
thus  stressing  the  importance  of  using  a carefully  designed  wave- launching 
structure. 

For  completeness,  the  solution  of  Eq.  (17)  illustrating  the  steady- 
state  mode-conversion  pattern  for  the  simple  case  of  a linear  density 
gradient  is  given  in  Appendix  A. 

Having  obtained  an  understanding  of  the  steady-state  response  of  the 

system  one  can  return  to  the  full  space-time  description  provided  by 

Eqs.  (14)  and  (15).  For  the  case  A = 0,  it  is  clear  that  these  equations 

decouple  into  two  independent  processes.  Eq.  (15)  describes  the  behavior 

of  local  (i.e.,  non-propagating)  electric  field  oscillations  in  the  direction 

perpendicular  to  BQ  . These  local  oscillations  are  not  short  circuited  in 

either  u>  < ui^^x)  or  oi  > Wpe(x).  Eq.  (14)  describes  the  propagation  of 

electric  field  oscillations  parallel  to  Bq.  These  oscillations  propagate 

up  to  the  oj  = u)pg(x)  point  and  are  reflected  back  to  the  outside  of  the 

plasma.  However,  for  A > 1 these  two  independent  oscillations  couple  to 

give  rise  to  a net  propagation  of  energy  beyond  the  u = u (x)  point. 

pe 

The  crucial  aspect  of  the  propagation  consists  of  the  excitation  of  a 
localized  oscillation  in  E^  by  the  gradient  of  E,.  Once  this  local 
oscillator  is  excited  at  a particular  point  in  space,  it  permits  the  pro- 
pagation of  subsequent  energy  pulses.  The  resulting  physical  picture  is 
that  of  a bootstrap  mechanism  whereby  the  Poynting  vector  of  the  incident 
radiation,  which  is  initially  aligned  with  Bq,  is  progressively  rotated 
through  a finite  angle  in  the  direction  perpendicular  to  Bq. 

The  relation  between  the  time  evolution  of  the  RF  energy  content  of 
the  system  and  the  Poynting  vector  can  be  found  from  Eqs.  (14)  and  (15). 


Multiplying  Eq.  (14)  by  E,*,  Eq.  (15)  and  E and  subtracting  the  complex 

Z X 

conjugate  of  the  resulting  expression  yields 


J 
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i 3t[|ExI2  * |EZI2J  * (E*3f  - c.c.) 


♦ 2i  [Im(e22)  JEz|2  ♦ hfej  |Ej2] 


(20) 


= 2iA  Re[3x(E*Ex)] 


Next , integrate  both  sides  of  Eq.(20)  from  an  initial  point  £ located 
outside  the  plasma,  up  to  a point  deep  inside  the  plasma.  At  ? , 

Ez  is  assumed  to  be  fixed  by  the  wave- launching  structure,  while  at 

Ex  “ E2  " 0 t»ecause  the  RF  energy  has  not  had  enough  time  to  propagate  to 

this  point.  Performing  the  integration  yields  an  expression  of  the  form 


i 


3_ 

3 1 


U + i 2k  U = I.  + I_ 

o a)  12 


(21) 


where. 


U = 


+ 


(22) 


represents  the  instantaneous  electric  field  energy  per  unit  area  stored  in 
the  region  bounded  by  the  points  xq  and  xf  (corresponding  to  ?q  and  £f). 

In  Eq.  (21) 


<y> 

w 


rx. 


|EXI2  * ImU„) 

0 


(23) 
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represents,  physically,  the  average  rate  of  energy  absorption  or  heating  which 
takes  place  between  xq  and  x^.  The  remaining  terms  in  Eq. (21)  are 


I 2 - 

I - - f—  dx  (E  * 3 E - C.C.) 
Ik  z x z 

o 

X 

o 


which  upon  integration  by  parts  yields 


I.  ■ -r—  (E  * 3 E - C.C.) 
1 k z x z 
o 


I,  * Im(E  * 3 E ) 
lk  z x z x 

o o 


similarly. 


I-  = iA  dx  (E  * 3 E + E * 3 E + C.C.) 
2 z x x x x z 


hence. 


I.  =*  -2iA  Re(E  E *) 

2 x z x 


Combining  Eqs.  (21)  through  (26)  one  obtains 

-rXRe(ExV>,  • (28) 

k oo  o 

o 

The  terms  on  the  left-hand  side  of  Eq.  (28)  can  be  readily  interpreted;  they 
simply  represent  the  rate  of  change  of  the  total  electric  field  energy  per 
unit  area  contained  in  the  system.  The  meaning  of  the  right-hand  side  of 
Eq.(28)  is  made  clear  by  evaluating  the  Poynting  vector  S at  xq 

4i  S =■  c Re  (E  x B*)  (29) 


B = V x E/(ik  ) 
~ ~ o 


hence 


I 
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B - ?(ikE  - 3 E )/(ik  ) 
J X x z o 


(31) 


where  y is  the  unit  vector  in  the  y direction.  Accordingly, 


thus. 


but. 


4rrS  = t—  Re[(i3  E * - kE  *)  (E  x - E £)] 

^ « X Z X Z X 

o 


4itS-x  = — [Re(iE  3 E *)  - k Re  (EE*)  ] 

K Z X Z Z X 

o 


(32) 


(33) 


Re(iE  3 E *)  = - Im(E  3 E *) 
z x z z x z 


(34) 


hence, 


4ttS  =■  - — Im(E  3 E *)  - CX  Re(E  E *) 
x k z x z ' zx 

o 

which  exactly  reproduces  the  right-hand  side  of  Eq.  (28)  when  evaluated  at 
x * x^.  The  interpretation  of  Eq.  (28)  is  now  complete,  it  reads  simply 


(35) 


•jr  U = 4tt(S  ) 
dt  xx 


(36) 


thus  stating  that  the  increase  in  the  electric  field  energy  content  of  the 
plasma  arises  due  to  the  development  of  a non-zero  just  outside  the  plasma 
edge. 

To  understand  the  nature  of  the  external  radiation  in  this  problem  it 

is  useful  to  consider  the  system  in  the  absence  of  plasma.  In  this  case, 

e = e ■ 1,  and  e * 1 - X2  <0.  Hence,  Eq.  (17)  takes  the  form 
zz  xx  ’ ^ 


3c-2E  - (X2-1)E  = 0 

c,  z z 


(37) 


whose  solution  is 


E » E exp(-  [X2-l]^“  k x) 
z o o 


(38) 


thus  indicating  that  the  external  radiation  Is  evanescent  in  the  direction  in 
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which  the  plasma  density  increases  (i.e.,  when  the  plasma  is  present).  In 
Eq.  (38)  Eq  is  a real  coefficient  which  represents  the  peak  field  produced  by 
the  wave-launching  structure  in  vacuum. 

The  energy  flow  associated  with  the  external  radiation  is  found  by  evalu- 
ating the  Poynting  vector  from  Eq.  (32),  in  which 


E - 3_E 

X 1-X2  ? 2 


hence. 


4 ITS  = — Re 
~ k 

o 


i 3 E 

— (xE  - zE  ) 

1-X2 


and  since  E^  is  purely  imaginary  and  given  by 


iX 


x [X2-l]1/2  x 


it  follows  that 


(39) 


(40) 


(41) 


4irS  = E 2 exp[-2(X2-l)1/2  k x]  . (42) 

2 X -1  ° ° 

Equation  (42)  illustrates  that  the  external  radiation  transports  energy  only  in 
the  direction  of  Bq  and  that  the  magnitude  of  this  flow  decreases  exponentially 
in  the  direction  of  increasing  plasma  density.  Accordingly,  the  picture  of 
the  system  at  a time  shortly  before  mode-conversion  occurs  is  as  illustrated 
qualitatively  in  Fig.  2. 

The  nature  of  the  oscillations  excited  in  the  mode- conversion  process  can 
be  better  understood  by  examining  the  dependence  of  the  generated  on  the 
phase  of  the  mode-converted  waves.  This  relationship  is  simpler  to  examine 
for  the  steady-state  case  in  which  the  electric  field  components  have  the 


form 
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Ez(x)  - a(x)  exp[ia(x) ] 


Ex(x)  * 6(x)  exp[ib(x) ] 


where  (a,  g)  and  (a,  b)  are  real  functions  describing  the  amplitude  and  phase 
of  the  (E  , E ) components.  Using  Eq.  (31)  leads  to 

Z X 


4ir  S - t — Re[i(3  a - ia3  a)  e - kg  e 
x k x x 

o 

c 2 

4tt  Sx  “ — [a  3 a - kag  cos(a-b)] 
o 

2 

Again,  in  the  region  of  interest  e = 1,  e - 1-A  hence 

XX  XX 

E - -i*-  i-  i-E 
x . ,2  k 3x  z 
1-A  o 

implying  that 


A I fl  ./  i -.2  .2.1/2 

7T~  ~ t(a3xa)  + (3x°°  ] 

A -1  o 


and  b-a  = <f>,  with 


4)  = tan 


-1  -3xa 


thus  Eq.  (44)  takes  the  form 


4*  sx  - {a23  a - cos<t>  [(c^a)2  + Oxct)2]1/2} 

o A -1 


using  Eq.  (47)  leads  to 


c 2 A 

4tt  S - -r — a 3 a 1 ~ — 

x k x ,2  , 

o A -1 


, c 2 9xa 

4’ s*  ’ - r “ XT 

o A -1 

Since  A2  > 1 to  satisfy  the  accessibility  criterion,  Eq.  (49)  in  conjunction 
with  Eq.  (36)  implies  that  in  order  to  couple  RF  energy  to  the  interior  of 
the  plasma  the  oscillations  excited  by  the  mode-conversion  process  must  be 
strictly  of  the  backward  type,  i.e., 
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3xa  < 0 . . (50) 

In  principle,  Eqs.  (14)  and  (15)  can  be  transformed  into  an  integral 
representation  by  using  a double  Laplace  transform  in  t and  £.  However,  the 
resulting  expressions  are  not  of  much  use  when  actually  trying  to  display  the 
space-time  evolution  of  the  electric  field.  In  particular,  this  task  is  rather 
arduous  for  complicated  density  profiles  and  is  only  limited  to  linear  problems. 
A more  direct  approach  is  to  solve  Eqs.  (14)  and  (15)  by  using  difference 
schemes  in  a digital  computer.  In  this  manner  one  can  handle  arbitrary  density 
profiles,  complicated  boundary  conditions,  and  in  addition  it  permits  the 
investigation  of  important  nonlinear  effects.  The  implementation  of  such  a 
numerical  study  for  Eqs.  (14)  and  (15)  is  discussed  briefly  in  Sec.  III.  For 
completeness,  the  integral  representation  of  Eqs.  (14)  and  (15)  is  given  in 
appendix  B for  the  simple  case  of  a constant  density  profile. 

The  space-time  formulation  provided  by  Eqs.  (14)  and  (15)  permit  the 

description  of  numerous  linear  and  nonlinear  processes  which  arise  due  to  the 

explicit  dependences  of  E and  e on  the  plasma  parameters.  Some  of  these 

XX  z z 

parameters  may  contain  interesting  zero-order  fluctuations  of  their  own,  or 
can  be  altered  in  a nonlinear  manner  when  the  applied  radiation  is  rather 
intense.  An  interesting  problem  associated  with  natural  fluctuations  is  the 
purely  linear  backscattering  of  the  evanescent  radiation  by  surface  ripples. 

This  is  a process  which  may  inhibit  the  rotation  of  S by  a significant  angle, 
however,  its  investigation  is  outside  the  scope  of  the  present  study. 

Out  of  the  many  nonlinear  processes  that  may  arise  at  the  edge  of  the 
plasma  when  large  power  levels  are  applied  in  a heating  experiment,  in  the 
present  study  we  have  chosen  to  isolate  the  specific  effects  associated  with 
the  modification  of  the  density  profile  by  the  ponderomotive  force.  This  is 
an  effect  whose  consequences  may  be  of  importance  because  it  alters  directly 
the  process  of  mode-conversion  by  shifting  the  point  where  ui  - Wpe(x)  and 
simultaneously  changes  the  gradient  scale  length.  A3  it  has  been  seen  earlier 
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in  this  section,  the  mode-conversion  process  is  responsible  for  the  generation 

of  a non-zero  Sx>  thus  the  profile  modification  can  play  an  inhibiting  role  in 

the  coupling  of  RF  energy  at  large  ppwer  levels.  It  is  also  of  general  interest 

to  compare  the  effects  of  the  nonlinear  density  changes  in  the  present  magnetized 

12  13 

plasma  problem  with  those  previously  observed  experimentally  and  theoretically 
in  the  interaction  of  intense  evanescent  radiation  in  an  unmagnetized  plasma. 

For  the  sake  of  simplicity,  in  the  present  study  we  introduce  the  nonlinear 
profile  modification  explicitly  through  the  expression  for  the  plasma  density,  i.e.. 


by  using 


n(x,  | E I ) = nQ(x)  exp  - 4—^ 

o 


in  which  nQ  represents  the  spatially  dependent  zero-order  density  (i.e.,  unmodified), 
T is  the  plasma  temperature,  and  |e|  is  to  be  determined  self-consistently  from 
Eqs.  (14)  and  (15).  The  nonlinear  density  shown  in  Eq.  (51)  appears  in  Eqs.  (14) 
and  (15)  through  the.  expression  = 4Tre  n(x,  |e|  )/m  contained  in  the  dielectrics. 

Clearly,  Eq.  (51)  is  an  oversimplification  of  the  actual  behavior  of  the 
plasma  density  which  is  subjected  to  the  intense  power  levels  expected  in  a fusion- 
grade  experiment.  Immediately,  in  a problem  of  this  nature  one  can  think  of  includ- 
ing a variety  of  effects,  such  as  the  generation  of  ion  acoustic  waves,  shocks, 
plasma  heating,  etc.  However,  before  proceeding  to  include  such  complications  one 
must  assess  what  results  arise  from  the  simple  model  provided  by  Eq.  (51).  Such 


is  the  philosophy  behind  the  nonlinear  studies  presented  in  Sec.  IV. 
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III.  NUMERICAL  PROCEDURE 

In  Implementing  the  numerical  study  of  the  mode-conversion  process  one 
proceeds  first  to  discretize  the  time  variable  in  finite  increments  At,  and  to 
time-average  both  Eqs.  (14)  and  (15)  to  obtain 

£ 


~ (F-?)  + j d^2(F+t)  + (F+?)  = 3 (G+tf) 

~ (G -4)  + | (G-Hf)  - ~ 3^(F +f) 


(52) 


(53) 


where  F - E2(£>  T) , G ■ Ex^>  T)  represent  the  present  values  of  the  electric 
field  components,  and  f,  # refer  to  their  value  at  the  previous  time  step,  e.g., 
4 * A(£,  x-Ax).  Eqs.  (52)  and  (53)  provide  a coupling  between  F and  G which 
must  be  solved  implicitly  at  every  time  step. 


Eliminating  G from  Eq.  (53)  gives 

«[(i/Ax)  - (e/2)]  + i | 3?(F+F) 
G " i/Ax  + e/2 


(54) 


which  can  be  inserted  in  Eq.  (52)  to  yield 


2 [(i/Ax)  + (e/2)] 

£ + [(i/Ax)  + (1/2)] 


2X 

Ax 


[(i/Ax)  + (e/2)] 
+ [(i/Ax)  + (1/2)] 
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Ax 


(55) 


in  which  one  makes  explicit  use  of  the  simplifying  properties  of  the  edge  region, 

namely  e^  =1.  Eq.  (55)  provides  a second  order  ordinary  differential  equation 

for  the  present  value  of  E^.  This  equation  can  be  solved  implicitly  if  one 

knows  the  previous  value  of  E^,  the  spatial  dependence  of  the  dielectric  e , 

and  the  appropriate  boundary  conditions. 

To  solve  Eq.  (55)  one  goes  through  an  averaging  and  discretization  in  £, 

similar  to  that  used  in  generating  Eqs.  (52)  and  (53).  The  resulting  second 

14 

difference  equation  in  £ is  solved  implicitly  by  using  the  well-known  tri- 
diagonal algorithm.  After  having  calculated  F by  this  method  one  returns  to 
Eq.  (54)  to  update  G,  thus  completing  one  full  time  step  in  the  problem. 
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The  boundary  conditions  used  in  solving  the  present  problem  are  the 
following:  1)  Et  is  held  fixed  at  the  value  Eq  at  a position  well  outside 

the  plasma,  and  2)  Ez  = 0 at  the  boundary  located  deep  inside  the  plasma. 

The  exact  location  of  the  point  outside  the  plasma  has  been  varied  for  a 
given  density  profile  to  verify  that  the  results  are  independent  of  its 
particular  numerical  value.  Similar  sensitivity  studies  have  been  made 
concerning  the  incremental  space-time  steps  A£  and  At . To  insure  that  no 
reflection  problems  develop  with  the  Ez  = 0 boundary  at  late  times  in  the 
evolution,  a strong  edge  absorption  is  included  through  e„,;  this  boundary 
absorption  does  not  give  rise  to  observable  effects  in  the  region  of  interest 
here. 

The  choice  of  boundary  condition  outside  the  plasma  allows  maximum 
freedom  to  the  mode-conversion  process,  since  the  value  of  3^Ez  is  left 
unspecified,  as  should  be  the  case.  This  derivative,  as  well  as  the  value 
of  Ex  at  the  outside  boundary  is  allowed  to  evolve  self-consistently  in  time, 
so  their  combination  gives  rise  to  the  gradual  development  of  a finite  S^. 

The  initialization  of  the  system  is  accomplished  by  assigning  the  steady- 
state,  pure  vacuum  solution  to  the  quantities  E„  and  E^ , i.e.,  they  are 
given  the  values  predicted  by  Eqs.  (38)  and  (39).  This  procedure  can  be 
interpreted  physically  as  a statement  that  prior  to  the  development  of  the 
mode-conversion  process  the  plasma  is  virtually  insensitive  to  the  RF  fields. 
However,  a few  cycles  after  turning-on  the  oscillator  the  plasma  begins  to 
respond,  and  it  is  this  time  that  we  arbitrarily  denote  as  t = 0.  Of  course, 
one  can  also  consider  the  slow  turn-on  problem  in  which  one  legislates  Eq  to 
be  a slowly  increasing  function  of  time,  however  no  additional  physical 
insight  is  gained  by  this  procedure  which  introduces  further  arbitrary  choices 
of  parameters. 

Various  shapes  of  the  zero-order  density  profile  nQ(x)  have  been 
investigated  in  which  in  general  the  density  is  taken  to  increase  monoton ically 
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as  a function  of  x.  Although  for  each  individual  profile  the  spatial  scale 
of  the  fields  developed  is  observed  to  vary  in  the  expected  manner,  the 
general  features  presented  in  Sec.  IV  are  insensitive  to  the  specific  choice 
of  profile. 

Finally,  the  nonlinear  density  profile  given  by  Eq.  (51)  is  introduced 
into  the  numerical  scheme,  at  a given  time  step,  by  averaging  the  quantity 
| E | over  the  two  previous  time  steps.  In  this  manner  a certain  amount  of 
numerical  inertia  is  retained  in  the  method. 

We  proceed  next  to  discuss  the  results  obtained  by  implementing  this 
numerical  scheme  with  the  help  of  the  mathematical  on-line  system  at  the 
University  of  California,  Los  Angeles. 
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IV.  RESULTS 

A.  Linear  Behavior 

In  the  investigation  of  Eqs.  (14)  and  (15)  a particular  zero-order  density 
profile  must  be  chosen.  In  general,  we  characterize  the  profiles  used  by  means 
of  a single  scale- length  L,  i.e.,  nQ  = nQ(x/L) . Thus,  it  proves  convenient 
for  display  as  well  as  for  physical  reasons  to  scale  the  spatial  variable  in 
terms  of  this  natural  length  of  the  problem.  This  choice  implies  that  in  the 
linear  regime  the  system  depends  explicitly  on  two  lumped  parameters,  namely 

kQL  and  A ==  k/kQ.  The  results  to  be  presented  correspond  to  the  particular 

2 

choice  A =2,  which  is  representative  of  the  expected  experimental  values. 

In  addition,  the  parameter  regime  of  interest  to  the  physics  of  the  plasma 
edge  corresponds  to  kQL  ~ 1 . 

To  give  a feeling  for  the  numbers  that  these  choices  may  imply  for  actual 
applications,  it  should  be  noted  that  for  incident  radiation  having  a frequency 
u = id/2tt  ' 3 GHZ,  which  may  be  considered  typical  of  a fusion-grade  experiment, 
the  length  scale  corresponds  to  L ' 1 cm.  However,  for  the  purposes  of  the 
present  study  no  specific  values  of  w and  L need  to  be  assigned. 

A characteristic  density  profile  used  in  the  present  studies  is  indicated 
by  the  solid  curve  shown  in  Fig.  3.  It  corresponds  to  the  choice 


n (x) 

-2 - [1  + tanh  (£  - 4)]  (x/2L) 


(56) 


where,  np  represents  the  density  for  which  Upg(x)  = u>,  and  whose  numerical 
value  does  not  need  to  be  specified.  The  dashed  curve  in  Fig.  3 is  the 
dielectric  e ; its  zero-crossing  separates  the  evanescent  and  propagating 
regions.  The  dotted  curve  in  Fig.  3 is  the  quantity  | E | 2 prior  to  the 
development  of  the  mode-conversion  process;  its  strong  evanescence  into  the 
plasma  is  apparent. 


W 
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It  should  be  emphasized  that  the  results  to  be  presented  in  the  remainder 
of  this  subsection  are  strictly  linear,  i,e.,  the  nonlinear  density  profile  given 
by  Eq.  (51)  is  replaced  by  nQ,  the  zero-order  unperturbed  profile. 

The  space-time  evolution  of  the  process  of  linear  mode-conversion  is 
illustrated  in  Fig.  4,  in  which  we  plot  the  spatial  dependence  of  the  real  and 

imaginary  parts  of  E^  for  various  different  times.  These  results  are  obtained 

2 

for  the  profile  shown  in  Fig.  3,  and  the  values  kQL  = 1,  (k/kQ)  =2.  It  should 
be  noted  that  there  is  a change  in  the  spatial  grid  between  the  displays  of 
Figs.  3 and  4.  In  Fig.  4(a)  ut  = 1.0,  hence  the  early  pattern  observed  is 
essentially  the  evanescent  vacuum  solution  except  for  the  appearance  of  the 
small  imaginary  part  near  the  edge.  By  ut  = 5.0  it  is  seen  in  Fig.  4(b)  that  the 
small  imaginary  part  of  Fig.  4(a)  has  grown  in  amplitude  and  in  doing  so  it 
becomes  spatially  out  of  phase  with  the  real  part.  The  combined  real  and 
imaginary  parts  have  developed  by  this  time  almost  a full  wavelength  which 
extends  into  the  plasma  beyond  the  point  where  e = 0.  In  Fig.  4(c)  one  clearly 
observes  the  development  in  the  interior  of  the  plasma  of  a wave  which  is 
continuously  linked  to  the  evanescent  external  radiation  in  the  region  x/L  ~ 3 
(i.e.,  near  zzz  = 0).  Fig.  4(d)  at  ut  = 20.0  indicates  the  existence  of  a 
well-defined  propagating  wavepacket  having  a significant  amplitude  and  whose 
role  is  to  transport  the  RF  energy  to  the  interior  of  the  plasma. 

The  transport  of  RF  energy  associated  with  the  sequence  of  events  illus- 
trated in  Fig.  4 can  be  seen  in  Fig.  5.  In  this  figure  we  reproduce  the 
information  of  Fig.  3,  but  include  in  it  the  spatial  dependence  of  the  mode- 
converted  quantity  [E_|2for  ut  = 20.0  (represented  by  the  dotted  curve).  For 
reference  it  should  be  mentioned  that  in  the  absence  of  mode-conversion  the 
quantity  jE,|2  would  be  evanescent  beyond  the  e,,  = 0 point. 

Although  the  sequence  of  events  depicted  in  Fig.  4 clearly  indicate  that 
RF  energy  is  being  carried  away  from  the  edge  toward  the  interior  of  the  plasma, 
it  is  not  clear  from  this  presentation  if  the  excited  mode  corresponds  to  a 

L J 
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forward  or  to  a backward  wave.  To  check  on  this  feature  one  needs  to  plot  a 
phasor  diagram  akin  to  that  commonly  used  in  sampling-time  interferometric 
measurements  of  plasma  waves.  Such  a plot  is  given  in  Fig.  6 for  the  E^ 
component  at  the  fixed  time  ut  = 5.0,  and  obtained  for  a profile  which  differs 
slightly  from  that  shown  in  Fig.  3.  In  Fig.  6 the  direction  of  the  arrows 
indicate  the  direction  of  increasing  density,  or  equivalently,  to  penetrating 
deeper  into  the  plasma.  The  fact  that  in  Fig.  6 the  phasor  rotates  in  the 
clockwise  direction  (i.e.,  in  the  direction  of  decreasing  phase  angle) 
identifies  the  type  of  modes  shown  in  Fig.  4 as  being  backward  waves,  as  is 
expected  from  the  analysis  leading  to  Eq.  (50).  The  corresponding  behavior 
for  the  Ex  component  is  illustrated  in  Fig.  7. 

It  has  been  demonstrated  in  Figs.  4 through  7 that  the  process  of  linear 
mode-conversion  indeed  takes  place  as  expected  and  that  it  generates  a backward 
wavepacket  that  proceeds  to  fill  the  interior  of  the  plasma  with  RF  energy. 
Having  observed  the  evolution  of  this  process  it  is  of  interest  to  inquire 
how  the  steady-state  is  attained  for  a given  spatial  point  in  the  region  where 
e„z  > 0.  The  answer  to  this  question  is  provided  by  the  sequence  of  curves 
shown  in  Fig.  8.  In  this  figure  we  plot  the  total  electric  field  energy  density 
|Ej2  as  a function  of  spatial  location  with  time  as  the  parameter  being  varied 
over  the  interval  0.5$  vt  s 10.0.  Again,  this  result  is  obtained  for  a profile 
slightly  different  from  that  of  Fig.  3,  and  with  kQL  = 1,  (k/kQ)2  = 2.  The 
salient  feature  uncovered  by  Fig.  8 is  that  the  RF  energy  propagates  away  from 
the  e = 0 region  via  a leading-edge  pulse  which  grows  in  amplitude  spatially 
as  it  penetrates  deeper  into  the  plasma.  At  a fixed  spatial  location,  however, 
nothing  occurs  in  time  until  the  foot  of  the  pulse  arrives.  After  its  arrival 
at  a given  point  the  local  energy  density  grows  in  time  until  the  peak  of  the 
leading-edge  pulse  passes  by.  Shortly  after  the  passage  of  the  propagating 
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pulse,  the  amplitude  of  the  local  energy  density  stops  growing  and  attains  its 
steady-state  value.  The  attainment  of  the  steady-state  is  reflected  in  Fig.  8 
by  the  accurate  overlapping  of  all  curves  in  the  trailing  edge  of  the  propa- 
gating pulse. 

B.  Nonlinear  Behavior 

We  proceed  next  to  describe  the  results  obtained  when  the  nonlinear  density 
modification  given  by  Eq.  (51)  is  included  in  the  space-time  evolution  of  the 
mode-conversion  process.  A characteristic  profile  distortion  produced  at  the 
edge  of  the  plasma  is  shown  in  Fig.  9.  In  this  figure  the  dashed  curve  corres- 
ponds to  the  zero-order  unperturbed  density,  while  the  continuous  curve  is  the 
nonlinearly  modified  profile  at  ut  = 10.0,  for  kQL  = 1,  and  (k/kQ)2  = 2. 

Three  qualitatively  different  types  of  distortions  can  be  identified  in 

Fig.  9,  namely:  1)  the  recession  of  the  point  <d  (x)  = w away  from  the  wave- 

pe 

launching  structure;  2)  the  local  steepening  of  the  profile  at  the  newly  created 

a (x)  = u point;  and  3)  the  generation  of  density  pockets  or  striations  beyond 
pe 

this  point.  All  of  these  distortions  work  in  the  direction  of  decreasing  the 

effectiveness  of  the  mode  conversion  process,  hence  the  nonlinearity  in  this 

problem  is  self-limiting  in  the  sense  that  the  distortion  leads  to  a reduction 

of  | E | 2 inside  the  plasma  instead  of  enhancing  it.  This  behavior  is  quite  the 

opposite  of  that  observed  in  the  mode-conversion  problem  in  an  unmagnetized 

12  13 

plasma.  As  is  well  known,  ’ in  the  unmagnetized  case  the  generation  of 
density  cavities  localizes  the  electric  field  and  leads  to  the  formation  of 
high  intensity  spikes. 

The  quenching  effect  produced  on  the  mode-converted  wave  by  the  modifica- 
tion of  the  profile  is  illustrated  in  Fig.  10  for  a study  made  with  a zero- 
order  density  of  the  form  n^n^  = x/L,  and  with  kQL  = 1,  (k/kQ)2  = 2.  Fig.  10 

is  a phasor  plot  of  E_  (of  the  type  discussed  earlier  in  connection  with  Fig.  6) 

z 

in  which  we  present,  on  the  same  scale,  the  patterns  obtained  for  a strictly 
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linear  case.  Fig.  10(a),  and  for  a strongly  nonlinear  run.  Fig.  10(b).  The 
nonlinear  result  shown  is  obtained  with  an  initial  external  excitation  level 
(i.e.,  Eq  in  Eq.  (38)  that  produces  an  initial  value  of  | E| 2 / (.■X^nQT)  = 7 at  the 
point  where  wpe(x)  = 

It  is  clearly  seen  in  Fig.  10  that  the  nonlinearity  reduces  the  amplitude  of 
the  mode  converted  wave,  thus  producing  a significant  quenching  in  the  amount 
of  RF  energy  that  penetrates  to  the  interior  of  the  plasma. 

The  salient  features  observed  in  Fig.  10(b)  are:  1)  the  increased 

evanescence  (straight-line  pattern)  due  to  the  recession  of  the  cd  = (x)  point 

pe 

to  the  interior;  2)  the  appearance  of  an  abrupt  phase  change  associated  with  the 
increased  steepness  of  the  profile  (manifested  through  the  rapid  transition  from 
straight-line  to  a curve) ; and  3)  an  overall  retardation  in  the  evolution  of  the 
mode-converted  wave  relative  to  the  linear  case.  A similar  quenching  occurs 
also  in  the  E^  component,  as  is  illustrated  in  Fig.  11  for  a case  kQL  = 2, 

(k/kQ)2  = 2,  and  in  which  the  chosen  Eq  yields  | e| 2/(4wnoT)  = 0.5  at  the  Wpe  point 
initially. 

A convenient  method  of  assessing  the  overall  effect  of  the  nonlinearity 
consists  of  evaluating  the  instantaneous  RF  energy  content  of  the  plasma  slab 
considered,  i.e.,  by  calculating  the  normalized  quantity 


x2 


u(t) 


dx| E| 2 

. X1 

E2 


(57) 


o o 

where,  xj  and  X2  refer  to  the  boundaries  of  the  slab.  In  Eq.  (57)  the  division 
by  Eq2  removes  the  trivial  dependence  on  the  external  driver  amplitude,  there- 
fore any  changes  obtained  in  u(t)  by  varying  Eq  are  due  to  explicit  nonlinear 
changes  in  the  profile.  Such  a dependence  is  illustrated  in  Fig.  12  for  an 
initial  profile  n^/n^  = x/L,  kQL  = 1,  and  (k/kQ)2  = 2.  The  ton  curve  corresponds  to 
a strictly  linear  case  (i.e.,  the  nonlinearity  is  shut-off),  while  the  other  curves 
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correspond  to  initial  values  of  |E  | 2/ (j+irn^T)  at  the  ojpe  point  of  0.4,  1.6,  and 
3.6  respectively. 

It  is  evident  from  Fig.  12  that  as  one  increases  the  external  power  in  the 
wave-launching  structure,  the  relative  amount  of  RF  energy  which  actually 
penetrates  into  the  plasma  decreases  steadily.  This  result  can  be  interpreted 
in  the  practical  sense  as  saying  that  the  RF-plasma  coupling  at  the  surface  is 
a decreasing  function  of  the  applied  power.  Physically,  of  course,  this  means 
that  the  nonlinear  quenching  of  the  mode-conversion  process  decreases  the  angle 
through  which  the  Poynting  vector  can  be  rotated. 

From  the  results  of  Fig.  12  one  can  extract  the  rate  at  which  the  plasma 
slab  is  filled  with  RF  energy,  i.e.,  by  calculating  the  quantity 


V = d_  u(t) 
c “ dt  c 


(58) 


in  which  u(t)  is  defined  by  Eq.  (57)  and  the  speed  of  light,  c,  proves  to  be  a 
convenient  scaling  factor  since  the  quantity  V denotes  an  effective  speed  of 
energy  propagation  into  the  plasma.  Clearly,  the  quantity  V/c  can  also  be 
interpreted  as  the  scaled  Poynting  vector  perpendicular  to  Bq. 

The  amplitude  dependence  of  V/c  is  illustrated  in  Fig.  13  for  the  results 
previously  shown  in  Fig.  12.  The  accentuated  squares  in  Fig.  13  denote  the 
numerical  results,  and  the  continuous  curve  is  an  excellent  eyeball-fit  given 
by  the  expression 


7 = (0.11)  exp  [-(1.3)  /p] 


(59) 


where  p = | E | 2/ (4irnQT)  at  the  uipe(x)  = to  point  prior  to  mode-conversion. 

It  is  of  interest  to  inquire  what  effect  the  nonlinearity  has  on  the 

behavior  of  the  propagating  leading -edge  pulse,  observed  in  the  linear  regime 

and  whose  evolution  is  shown  in  Fig.  8.  A characteristic  nonlinear  response 

is  exhibited  in  Fig.  14,  in  which  we  present  the  spatial  dependence  of  | E | 2 and 

n for  different  times  in  the  interval  5.0  <vt  < 10.0,  for  the  case  k L = 1 

- - o’ 
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(k/kQ)2  = 2,  and  for  the  initial  profile  indicated  by  the  dashed  curve. 

It  is  seen  in  Fig.  14  that  the  leading-edge  pulse,  characteristic  of  the 
linear  regime , breaks  up  into  an  isolated  peak  which  proceeds  to  propagate  into 
the  plasma.  The  isolated  peak  observed  in  Fig.  14(c)  is  reminescent  of  a 
Langmuir-wave  soliton  and  its  associated  density  cavity.  However,  the  reason 
for  the  localization  of  the  pulse  in  Fig.  14(c)  is  entirely  different  from  the 
Langmuir-wave  case.  In  the  Langmuir-wave  soliton  the  localization  is  associated 
with  the  trapping  of  the  radiation  in  a density  cavity,  while  in  here  no  trapping 
phenomena  occurs.  In  here,  the  isolated  pulse  arises  because  its  leading  edge 
modifies  the  profile  thus  reducing  the  efficiency  of  mode-conversion  for  pulses 
that  trail  in  time.  In  other  words,  the  leading  edge  samples  the  higher 
efficiency  zero-order  profile,  while  the  trailing  edge  samples  the  modified 
profile  which  has  lower  efficiency  for  mode  conversion. 

The  nonlinearly  generated  RF  pulses  can  become  quite  dramatic  in  their 

appearance  for  shallow  density  profiles,  i.e.,  when  u>  (x)  does  not  become  too 

pe 

large  compared  to  w An  example  of  the  sharp  pulses  that  can  appear  in  these 
cases  is  shown  in  Fig.  15,  which  has  been  obtained  with  nQ/n  = 2 [l+tanh(£-4) ] 
for  kQL  = l,(k/k&)2  = 2,  and  |E | 2/ (dirn^T)  = 6.0  at  the  w = oipe  point,  initially. 
The  time  of  observation  in  Fig.  15  corresponds  to  ut  = 13.0. 

It  is  clear  from  the  nonlinear  physics  that  underlies  the  generation  of 
the  pulses  exemplified  in  Figs.  14  and  15,  that  multiple  peaks  or  pulse-trains 
should  appear  in  the  long  time  evolution  of  the  system.  The  reason  is  simply 
that  the  nonlinearity  is  self-limiting,  hence  after  the  passage  of  a large 
amplitude  peak  the  system  returns  to  its  linear  behavior  which  proceeds  to 
generate  another  linear  leading-edge  pulse,  and  thus  the  cycle  repeats.  The 
generation  of  multiple-peak  pulses  can  also  be  interpreted  by  saying  that  the 


x component  of  the  Povnting  vector,  S^,  is  constantly  oscillating  between  its 
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linear  and  nonlinear  values,  thus  leading  to  an  intermittent  flow  of  RF  energy 
into  the  plasma. 

The  existence  of  these  multiple-peak  pulses  or  pulse  trains  is  indicated  in 

Fig.  16,  in  which  we  display  the  spatial  dependence  of  | J 2 for  various  times 

in  the  interval  10.0  < ut  < 14.0  and  obtained  for  a case  n /n  = x/L,  k L = 1, 

- o p o 

(k/kQ)2  = 2,  and  jE | 2/ ^irn^T)  = 3.0  at  m = a>pe(x),  initially.  In  Fig.  16(a)  one 
observes  a dominant  single  peak  located  near  the  left  boundary  and  followed  by 
two  small  pulses.  Fig.  16(b)  shows  that  the  pattern  of  Fig.  16(a)  propagates 
into  the  plasma,  with  the  trailing  pulses  becoming  more  visible.  Finally,  in 
Fig.  16(c)  one  observes  a fully  developed  multiple-peak  pulse  train  that 
transports  the  RF  energy  to  the  interior  of  the  plasma. 
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V.  SUMMARY 

In  the  present  investigation  the  use  of  a modulational  representation  for 
the  electric  field  has  permitted  the  full  description,  in  the  linear  as  well  as 
in  the  nonlinear  regime,  of  the  process  of  mode-conversion  encountered  at  the 
edge  of  the  plasma  in  the  lower-hybrid  heating  scheme. 

It  has  been  found  that  the  salient  feature  of  the  mode-conversion  process 
consists  of  the  gradual  excitation  of  a backward  wave  which  causes  a rotation  of 
the  Poynting  vector  of  the  external  evanescent  radiation  in  the  direction  of  the 
density  gradient,  thus  allowing  the  RF  energy  to  penetrate  into  the  plasma. 

The  transient  transport  of  RF  energy  away  from  the  plasma  edge  toward  the 
interior  has  been  found  to  be  accomplished  by  the  propagation  of  a leading-edge 
pulse  whose  arrival  at  a given  spatial  location  denotes  the  onset  of  the  steady- 
state. 

It  has  been  found  that  the  nonlinear  profile  modification  produced  by  a 
simplistic  mode  of  the  ponderomotive  force  can  cause  a significant  quenching  of 
the  mode-conversion  process,  thus  reducing  the  amount  of  RF  energy  that  can  be 
coupled  to  the  interior  of  the  plasma  at  large  power  levels.  These  results  thus 
warrant  future  studies  in  which  more  complicated  and  realistic  models  of  the 
profile  modification  are  used. 

In  the  strongly  nonlinear  regime  the  present  study  has  uncovered  new  and 
interesting  multiple-peak  pulse  trains  that  can  be  quite  striking  in  shallow 
density  profiles.  The  consequence  of  the  existence  of  these  nonlinear  pulses 
is  that  the  flow  of  RF  energy  to  the  interior  of  the  plasma  is  intermittent  at 
large  power  levels. 


L 
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Accordingly,  the  steady-state  linear  mode-conversion  has  the  solution 


where  k - (X2-!)1^3  (k  L)2^3,  and  the  remaining  constant  y is  determined  by 


the  amplitude  of  the  external  radiation. 
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APPENDIX  B 

Defining  A = E^,  B = Eqs.  (14)  and  (15)  become 

n A + 9_2  A + e A-  iX3-B 
T C ZZ  £ 

i 3t  B + e B = iX3^A  . (B-l) 

Introducing  the  Laplace  transform  operator  in  time 

00 

s =^(s;t)  { } ■ | dT  eST  { } (B-2) 

0 

and  A = A,  B = B,  yields 
s s 


e+is  1 

(e  +is)A  = i 
zz 

e+is 

e +is 

e +is 

XX  J 

XX 

/ 

A(£,t=0) 


iX 


e +is 


XX 


3?B(C,T=0) 


(B-3) 


Further  introduction  of  a Laplace  transformation  in  space,  in  the  manner 
shown  in  (B-2),  i.e.,^^  = <*^(k;£)  and  = A transforms  (B-3)  to  the 

form 


[k2+K2(s)]A 


— X 

|3£  A 


(£=0,s)  + kA(C=0,s)  + i 


e+is 


e +is 
xx 


A(k, x=0)  - 


XX 


% 3EB<k'r-0) 


(B-4) 

2 

where,  K (s)  = (e+is)  (ezz+is) / (e^+is)  and  it  is  explicitly  assumed  that  the 
dielectrics  are  spatially  independent.  If  a spatial  dependence  were  allowed, 
then  (B-4)  would  become  an  integral  equation  involving  appropriate  convolutions 
of  the  dielectrics  and  fields. 

Using  the  pure  vacuum  radiation  solutions  to  evaluate  A(£,t=0)  and 
B(£,t=0)  through 


A =*  a exp[-(  |e/£xx|)1,/2  £] 


B 


( I e ! 


iX 


A 


(B-5) 


where  a denotes  the  peak  field  amplitude  produced  by  the  wave-launching 
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structure,  yields 


A(x,t) 


i ds 

r j,  st  kx 

. dk  e e r 

ri_  rl 

j 2iTi  ; 

2iTi  ,2  ^ 2 . . * 

k + K (s) 

t9C  SJ 

(£“0,s)  + kA(£>0,s) 


+ 


a 

(exx+is) (k  + [ *e/exxl ]1/2) 


(e+is) 


} 


(B.  6) 


where  j refers  to  the  appropriate  Laplace  inversion  contours. 

Equation  (B-6)  gives,  in  principle,  the  space-time  evolution  of  the 
mode-conversion  process  provided  the  value  of  A and  its  spatial  derivative  are 
known  at  the  edge.  The  simplest  response  contained  in  (B-6)  is  that  due  to  a 
fixed  amplitude  at  the  boundary.  Denoting  this  part  of  the  solution  as  A^ 
yields 


AIX(x,t) 


<a  , ds 
' 2 s 


St 


jgiK(s)x  + e-iK(s)Xj 


(B-7) 


from  which  one  obtains  the  asymptotic  solutions  (t  ®)  corresponding  to  the 
mode-converted  waves  of  the  form 


exp{±i  f I £ezz/£xx  1 3 1/2  5) 


(B-8) 
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FIGURE  CAPTIONS 

FIG.  I Cartessian  coordinate  system  and  corresponding  slab  geometry  used 
in  the  problem. 

FIG.  2 Qualitative  sketch  indicating  the  direction  of  power  flow  associated 
with  the  evanescent  radiation  prior  to  mode-conversion. 

FIG.  3 Spatial  dependence  of  typical  zero-order  quantities  investigated. 

The  continuous  curve  represents  the  density  profile,  the  dashed  curve 
denotes  the  dielectric  eT,,  and  the  dotted  curve  corresponds  to  | E | 2 
prior  to  mode-conversion. 

FIG.  4 Spatial  dependence  of  the  real  (continuous  curve)  and  imaginary 

(dashed  curve)  parts  of  E^  for  different  times  during  the  process  of 
mode- conversion.  Time  interval  displayed:  a)  ut  = 1.0,  b)  ut  = 5.0, 

c)  ut  = 10.0,  d)  ut  = 20.0. 

FIG.  5 Spatial  penetration  of  RF  energy  into  the  plasma  after  mode-conversion 
takes  place.  Continuous  curve  is  the  density,  dashed  curve  is  the 
dielectric  e^z-  The  dotted  curve  is  |Ej2  at  ut  = 20.0,  and  should 
be  compared  with  Fig.  3. 

/ 

FIG.  6 Phasor  plot  of  E^  at  ut  = 5.0.  Parameter  along  curve  is  the  spatial 

position  into  the  plasma,  and  increases  in  the  direction  of  the  arrows 
Clockwise  rotation  shows  mode-converted  wave  is  a backward  wave. 

FIG.  7 Phasor  plot  of  E^  at  ut  = 5.0.  Parameter  along  curve  is  the  spatial 

position  into  the  plasma,  and  increases  in  the  direction  of  the  arrows 
Clockwise  rotation  shows  mode-converted  wave  is  a backward  wave. 
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FIG. 


FIG. 


FIG. 


FIG. 


FIG. 


FIG. 


8 Space-time  penetration  of  RF  energy  in  the  linear  regime  over  the 
time  interval  0.5  < ut  < 10.0.  Overlapping  of  curves  for  different 
times  indicates  approach  to  the  steady-state. 

9 Typical  nonlinear  distortion  of  the  density  profile  obtained  for 
large  external  power  levels.  The  continuous  curve  is  the  profile 
at  ut  = 10.0,  and  the  dashed  curve  is  the  unperturbed  profile. 

10  Nonlinear  quenching  of  the  mode-conversion  process  observed  at  large 

power  levels.  Shown  aTe  phasor  plots  of  E7  for:  a)  linear  case,  and 

b)  strong  nonlinearity.  Curves  a)  and  b)  are  plotted  on  the  same 
scale  and  correspond  to  equal  times  in  the  evolution. 

11  Nonlinear  quenching  of  the  mode-conversion  process  observed  at  large 

power  levels.  Shown  are  phasor  plots  of  E^  for:  a)  linear  case,  and 

b)  mild  nonlinearity  (compared  to  that  used  in  Fig.  10).  Curves 

a)  and  b)  are  on  the  same  scale  and  correspond  to  equal  times  in 
the  evolution. 


12  Time  dependence  of  total  electric  field  energy  content  in  the  plasma, 
U,  for  different  values  of  the  externally  applied  amplitude,  Eq,  in 
arbitrary  units:  linear  case,  Eq  = 1.0,  Eq  = 2.0,  Eq  = 3.0. 
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Dependence  of  the  rate  of  electric  field  energy  penetration  into  the 
plasma 


1 

<1 

III 

0-1  CL 
*1 

x2 

dx|E/Ej2 

w l 

JX1  ' 

on  the  externally  applied  power  level,  in  arbitrary  units,  c is  the 
speed  of  light.  Heavy  squares  are  numerical  results,  and  solid 
curve  is  eyeball  fit  given  by  Eq.  (58). 
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FIG.  14  Nonlinear  space-time  penetration  of  RF  energy  into  the  plasma 
over  the  interval  5.0  < ut  < 10.0.  Dashed  curved  is  the  zero- 
order  density  profile;  continuous  curve  next  to  dashed  curve  is 
the  nonlinearly  modified  profile.  To  be  compared  with  Fig.  8 


FIG.  15 


FIG.  16 


Spectacular  nonlinear  pulses  and  density  striations  observed  at 
large  power  levels  in  shallow  density  profiles. 

Nonlinear  propagation  of  multiple-peak  pulse  trains  that  give  rise 
to  intermittent  transport  of  RF  energy  to  the  interior  of  the 
plasma.  Behavior  to  be  contrasted  with  linear  result  shown  in 
Fig.  8. 
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